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\ We analyze the nMSSM with CP violation in the singlet sector. We study the 

o 

static and dynamical properties of the electroweak phase transition. We conclude 

that electroweak baryogenesis in this model is generic in the sense that if the present 

limits on the mass spectrum are applied, no severe additional tuning is required 
Ph. 

to obtain a strong first-order phase transition and to generate a sufficient baryon 



asymmetry. For this we determine the shape of the nucleating bubbles, including 
the profiles of CP-violating phases. The baryon asymmetry is calculated using the 
advanced transport theory to first and second order in gradient expansion presented 
recently. Still, first and second generation sfermions must be heavy to avoid large 
electric dipole moments. 

PACS numbers: 98.80.Cq, ll.30.Er, 11.30.Fs 



Stephan. Huber @cern . ch 



Konstand@kth.se; 



T.Prokopec@phys.uu.nl 



M . G . Schmidt @ thphys . uni-heidelberg. de 



2 



I. INTRODUCTION 

Lately, electroweak baryogenesis (EWBG) [1] has again attracted more attention, not 
least because new collider data will hopefully provide information about the relevant (su- 
persymmetric?) physical parameters in the next years. 

Electroweak baryogenesis relies on a strong first-order electroweak phase transition as the 
source of out-of-equilibrium effects. During a first-order phase transition, bubbles of the low- 
temperature (broken) phase nucleate and expand to fill all space. An important aspect in the 
determination of the baryon asymmetry is the impact of transport phenomena . Without 
transport, CP- violating currents would only be generated near the bubble wall profile of 
the Higgs vevs. Close to the wall, sphaleron transitions are already strongly suppressed 
by the mass of the W-bosons, so that B-violating processes are inefficient in producing the 
observed baryon asymmetry (BAU). Early works that incorporated transport effects into the 
EWBG calculus were based on the WKB approach jjj. In this framework, a CP- violating 
shift in the dispersion relation induces a force of second order in the gradient expansion 
in the Boltzmann equation and leads to CP-violating fermion densities in the symmetric 
phase. Later on, this formalism was applied to the MSSM, where CP violation results from 
mixing effects in the chargino sector |4|. In this context, the formalism had to be extended 
to the case of mixing fermions. In the MSSM, the second-order source is too weak to yield 
a successful baryogenesis 

One disadvantage of the WKB method is the neglect of dynamical flavor mixing effects. 
While the shift in the dispersion relation is due to mixing of left- and right-handed compo- 
nents of the fermions and already present in the one-flavor case, flavor mixing contributions 
have been completely neglected after a flavor basis transformation to the mass eigenbasis. 
A series of papers [7|, |8j aimed at improving on this point by including flavor mixing by 
using a perturbative expansion of the Kadanoff-Baym equations. Here, the deviations of 
the Green function have been interpreted as sources in the diffusion equation. This approach 
- like the WKB-approach — has the weakness that the transport equations only describe 
the dynamics of two classical quasi-particles in the chargino sector. CP violation is commu- 



nicated 
of Refs. 



rom the charginos to the SM particles by their interactions. Therefore, the authors 



7. 



used the Winos and Higgsinos as quasi-particles in the interaction basis, where 



the interactions take a particularly simple form. In the WKB-approach the natural choice is 
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the mass eigenbasis. This dependence on a flavor basis is unsatisfactory, especially since the 
flavor mixing CP-violating source vanishes in the mass eigenbasis completely. A numerical 
analysis, making use of the first-order source in the interaction basis in Refs. 0, Q] , leads to 
successful electroweak baryogenesis for a certain range in the parameter space, even though 
at least some fine-tuning is required to fulfill the electron electric dipole moment (EDM) 
constraints. 

Recently, some of the authors have derived semiclassical transport equations for the 
chargino sector from first principles . The derivation is based on the Kadanoff-Baym 
equations and does not depend on classical reasoning to fill the gap between CP violation 
and transport effects. Technically, the two main improvements on the resulting transport 
equations are independence of the flavor basis and the absence of the source strength ambi- 
guities. First, since the semiclassical transport equations describe the dynamics of a 2 x 2 
matrix in flavor space, the transformation properties of the transport equations under flavor 
basis changes are explicit, and no restriction to quasi-particles has to be used. Secondly, 
because the CP-violating sources appear naturally and uniquely as higher order terms in 
the gradient expansion of the mass background fields, there are no ambiguities in the source 
strength. The mixing effects lead to an additional force of first order of the gradient expan- 
sion in the transport equations. In contrast, in the work reported on in [8|, for dimensional 
reasons, the sources had to be multiplied by a typical thermalization time r, while in j?| the 
removal of the ambiguity in the source strength was based on the (classical) Fick's law. In 
the case of second-order effects, the first principle derivation confirms the WKB approach if 

nn 

the latter is handled carefully [101. 

Applying this advanced transport theory to the MSSM [6|, the quantitative analysis shows 
two distinct features, which are less definitive in the results of the former approach [?], 
First, mixing effects are strongly suppressed away from mass degeneracy in the chargino 
sector, |"T.|± — m ^±l > (20 GeV) 2 . So mixing in that sector is only effective if the a priori 
unrelated Wino mass parameter M2 and the Higgsino mass parameter /1 are tied together. 
Secondly, the produced BAU suffers from an exponential Boltzmann suppression in the case 
of heavy charginos. In this formalism, successful electroweak baryogenesis requires rather 
large CP violation in the chargino sector, sin(^c*p) > 0.25, even for the most favorable 



choice of the other model parameters. In comparison, the approach followed in Refs. 



7: 



leads to viable baryogenesis for less constrained chargino masses and CP violation of order 
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sm(5cp) > 0.1. Hence, if the advanced transport theory is used, not only the parameter 
space of viable baryogenesis is much more restricted, but also the maximally achievable 
BAU is smaller. Because of the necessity of large CP- violating phases, additional arguments 
(cancellations, or a large value of the CP-odd Higgs mass parameter) are required to suppress 
the electron EDM by a factor 5-6. On top of this, a light right-handed stop and a light Higgs 
are needed to allow for a strong first-order phase transition 12j (light stop scenario). Thus, 
electroweak baryogenesis in the MSSM is severely constrained. 

In this paper we study electroweak baryogenesis in a singlet extension of the MSSM, where 



the divergences of the singlet tadpole are tamed by a discrete R-symmetry [13|, [14], Jj], Il6 |. 
The R-symmetry is violated by the supersymmetry breaking terms. A singlet tadpole is 
then induced at some high loop order, which is too small to destabilize the weak scale, but 
large enough to evade the cosmo logical domain wall problem [171 ] . 

Our analysis supports the result of Ref. 18j that a strong first-order phase transition is 
quite generic, once experimental constraints on the Higgs and sparticle spectrum are taken 
into account. The phase transition is induced by tree-level terms in the Higgs potential 
without the need of a light stop. Going beyond Ref. fis| . we actually compute the baryon 
asymmetry and the bubble wall properties. We find that the observed baryon asymmetry 
can be produced with mild tuning of the model parameters. The first and second generation 
squarks and sleptons have to be heavy (a few TeV) to suppress the one-loop contributions 
to the EDMs. In contrast to the MSSM, there are no strong constraints from the two-loop 



EDMs, since tan {3 is usually small [19 



Several variants of MSSM singlet extensions have been studied in the 
respect to their impact on electroweak baryogenesis 



21 



22 



23 



24 



25 



26 



iterature with 



detailed of these studies is Ref. 



271 ] . The most 



261 ]. where a general singlet model without discrete sym- 



metries was considered. This general model supports electroweak baryogenesis in a large 
part of its parameter space. In the current work, the R-symmetry forbids a self coupling 
of the singlet, leading to a quite constrained Higgs and neutralino phenomenology. Still, it 
is encouraging to see that even this restricted framework allows for successful electroweak 
baryogenesis. 

The paper is organized as follows. In section [IT] we will present the model and clarify 
notation. In sections [IV] and |V] we will discuss the mechanism that drives baryogenesis in 
the nMSSM and the dynamics of the phase transition. In section I VI I numerical results will 
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be presented before we conclude in section IVHI 



II. THE NEARLY MINIMAL SUPERSYMMETRIC STANDARD MODEL 



A. The Model 



The notation in this section follows Refs. [15|, 118J, including, however, a generalization 
to an additional CP-violating phase in the singlet sector. The superpotential, including the 
multi-loop generated tadpole term is 

2 

WnMssM = ^SHi ■ H 2 - -j±S + y t Q ■ H 2 U C + ■ ■ • , (1) 

where the dots denote the remaining terms in the MSSM superpotential, Hi = (H^,H^), 
H 2 = (Hf, H$), S is the singlet superfield and, A ■ B = e ab A a B b = A 1 B 2 - A 2 B X . 
The tree-level potential consists of 

V = V F + V D + VU, (2) 
where, restricting to third generation quarks, 

V F = |Aif x • H 2 - ^f\ 2 + IXSH? + y t t L t* R \ 2 

+\XSH^ + y t b L t* R \ 2 + \\S\ 2 H\H 2 

+ \y t t* R \ 2 HlH 2 + \y t Q-H 2 \ 2 , 
Vd = °j(HtH 2 - H\Hi) 2 + ^W 2 | 2 , (3) 
V soit = 771^1^ + m 2 2 HlH 2 + m 2 \S\ 2 

+ (t s S + h.c.) + (a x SHi ■ H 2 + h.c.) 

+m 2 Q Q j Q + mlrltnl 2 + (a t Q ■ H 2 t* R + h.c). 

The Higgs sector of this potential has only one physical CP-violating phase, which after 
some redefinition of the fields can be attributed to the parameter t s . We assume that this 
phase is the only source of CP violation in the model, i.e. the gaugino masses and squark 
and slepton soft terms are taken to be real as well as the parameter A. 
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In the case when the squarks have vanishing vevs, the tree-level Higgs potential is 

V = m\E[E x + m\H\E 2 + m 2 s \S\ 2 + A 2 |#i • H 2 \ 2 

+X 2 \S\ 2 (HlH 1 + HlH 2 ) + 9 ^{H\H 2 - HlH,) 2 + 9 -^\HlH 2 \ 2 

o Z 

+ t s (S + h.c.) + {a x SH x ■ H 2 + h.c.) - m\ 2 (H x ■ H 2 + h.c). (4) 

We define the vevs as (Hf) = (/^e 291 , (H®) = 4> 2 e' tg2 , (S) = <f) s e iga . We use the gauge freedom 
to set (-^r) = an d choose the phase convention qi = q 2 = q/2. The absence of a charged 
condensate (H 2 ) = will be ensured by the positivity of the squared charged Higgs mass 
j3| that is determined in the numerical analysis. 
Using the definition of the (3 angle 

(j) 1 = cj) cos(/3), 02 = sin(/3), (5) 

and t s = \t s \ e lQt , the tree- level potential reads finally 

V = M 2 <f) 2 + m 2 s <J) 2 + 2 \t s \<f) s cos(g t + q 8 ) + 2 50 2 s + A 2 </> 2 2 + A 2 4 , (6) 

where we have used 

M 2 = m\ cos 2 (3 + m\ sin 2 (3 — m\ 2 sin 2(3 cos q, 
a = y sm2f3cos{q + q s ), 

A 2 = ^ sin 2 2/3 + % cos 2 2/5, (7) 
4 8 

to shorten the notation. 



B. Effective Potential at Zero Temperature 

At zero temperature we take into account, in addition to the tree-level potential, the 
Coleman- Weinberg one-loop contributions 



AV = 




where the two sums run over the bosons and the fermions with the degrees of freedom 
and gf respectively, and 

(9) 



m 



Km 2 ) = -r 
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We choose the renormalization point to be Q = 150 GeV in the DR-scheme and suitable 
counter-terms, such that the one-loop contributions to the potential preserve the location 
of the tree-level minimum. This leads to the following shifts in the bare mass parameters 

. 2 1 dAV . 2 1 dAV . 2 1 dAV . . 

20i d(j)i 202 902 20 s d(j) s 

As degrees of freedom of the relevant one-loop contributions, we take 

9w = 6, g z = 3, g t = 12, g h = g- tR = 6. (11) 

Contributions from the charginos and neutralinos are not taken into account. The masses 
used in the one-loop potential are listed in Appendix |A] The neutral Higgs masses are 
computed from the second derivatives of the one-loop potential. 



C. Effective Potential at Finite Temperature 



Taking into account temperature effects, we correct the effective potential by the thermal 
one-loop contributions, which read 



AV 1 



rp4 

2^2 



5>J + (mg/T 2 ) -J>J_K/T 2 

. b f 



(12) 



with the definitions 



POO 

■J±(y 2 ) = / dxx 2 log(l =f exp(— \/ x 2 + y 2 ))- 
Jo 



(13) 



Discussing a strong first-order phase transition in the MSSM, it is important to modify 
this expression by undergoing a two-step procedure, first deriving a 3D effective action and 
then treating this further with two-loop perturbation theory [28] or, more safely, with lattice 
numerical methods [291 ] . A simplified prescription, the above one-loop expression modified 
by "daisy" resummation, follows the same direction. It allowed the formulation of the 
postulate of a "light" stop, in order to obtain a strong phase transition, in a transparent 
way [l2j. In the nMSSM model the strong first-order phase transition should be triggered by 
the tree-level Lagrangian and we do not need this refined analysis. Just adding the "daisy" 
correction does not necessarily improve the analysis. 



8 



species 


mass bound 


charginos 


> 104 GeV 


neutralinos x° 


> 25 GeV 


charged Higgses 


> 90 GeV 


neutral Higgses H° 


> 114 GeV 



TABLE I: Mass constraints on the spectrum 
III. CONSTRAINTS AT ZERO TEMPERATURE 



Before we analyze the phase transition and compute the produced baryon asymmetry, 
we confront the model with constraints coming from collider physics. The present limits are 
summarized in Tab. HI 

In models with extended Higgs sectors, the Higgs couplings deviate from the SM values. 
Of particular importance is the ZZTii vertex, where Tii denote the neutral Higgs mass 
eigenstates, as computed from the one-loop potential. The size of this vertex is reduced by 
a factor 

& = cos(/3)0 a + sin(/3)0 i2 , (14) 

with respect to the SM value. The matrix Oik relates the mass eigenstates with the two 
CP-even flavor eigenstates S\,S2- In the CP- violating case, O is a 5 x 5 matrix without 
block-diagonal structure, and the special form of Eq. (I14p is due to our convention for CP- 
even Higgs states (see Appendix [A]). If the neutral Higgs mass is below the value given in 



Tab. [H the LEP bound translates into an upper bound on £, as given in Ref. 30] . 



We do not implement any constraints on the squark spectrum, but choose the following 



stop mass parameters as used in Ref. [181 ]: 



m l Q = mfj = 500 GeV, a t = 100 GeV. (15) 

The nMSSM suffers from a light singlino state. Because of the missing singlet self cou- 
pling, this state acquires its mass only by mixing with the Higgsinos. This is an important 
difee.ce tan m0r e g ene Ial singl e t m ode ls , such a S the oee — ie Re, Q. If this 
lightest neutralino has a mass m^o < m z /2, it contributes to the invisible Z width, leading 



to the constraint [3 



BR(Z - xY) 



9 2 (\Ui^-\U^\r m z 
An 24 cos 2 ( 



1 



1 w , 



3 

4m| X " 
mi 



< 0.8 x 10" 



(16) 



Here, U denotes the unitary matrix that diagonalizes the neutralinos as defined in Ap- 
pendix |A] and Tz = 2.5 GeV denotes the Z width. 

A light neutralino can be avoided if the Higgs singlet coupling A is taken to be large. For 
large values of A a Landau pole is encountered below the GUT scale. Avoiding this Landau 



pole requires tan/5 > 1.3 and A < 0.8 



181 ]. but we will also consider larger values of A. This 



can be motivated by the so-ca 
at some intermediate scale 
In Ref. 



31 



led "Fat Higgs" models, where the Higgs becomes composite 



32]. 



181 ] the model was further constrained by the relic neutralino density. Here 



we will not impose this constraint directly. However, from Ref. 18] we take the bound 
m^o > 25 GeV, which ensures that the dark matter density remains below the observed 
value. For much larger masses the relic neutralino density is quite small, so that neutralinos 
will only provide a fraction of the total dark matter in the Universe. 



IV. ELECTROWEAK BARYOGENESIS 

A. Sources of CP violation in the nMSSM 

As discussed earlier, in the nMSSM there is the possibility of additional CP violation in 
the singlet sector. As in the MSSM, the relevant source in the transport equations usually 
comes from the charginos, even though the neutralinos can in certain cases contribute sizable 
effects as well 33J. In the interaction basis, where the Higgsinos and Winos are the quasi- 
particles, their mass matrix takes the following form: 
















w 


° i 



M 2 

#i(z)e -i ^ 



i 2 (z)e 



H(z) 



(17) 



where z denotes the direction along which a nearly planar bubble wall is moving. Unlike 
what happens in the MSSM, the effective \i term acquires a z-dependence 



fx(z) = -\<p s (z) e^*). 



(18) 
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The leading contribution (to the left-handed current) to second order in the gradient 
expansion is proportional to (see Eq. (92) in Ref. 9], and Ref. js|) 

5 (2) ~ {m^'m - m^m", d k J cq } D , (19) 

where g eq denotes the zero component of the vector part of the chargino Green function 
in thermal equilibrium and in the interaction basis. The superscript D indicates that the 
diagonal entries in the mass eigenbasis are projected out following the conventions of Ref. 
The first-order sources that are used to calculate the BAU in the MSSM and nMSSM in 
Refs. jg] and become prominent when chargino mass eigenstates are nearly degenerate. 
Here these sources are generally expected to be suppressed with respect to the second-order 
source (119]) . because we are mainly interested in the generic non-degenerate case. Note that 
gradient expansion applies when the typical momentum of the particles is large with respect 
to the inverse wall thickness, k 3> i' 1 . Since k ~ T, this condition is reasonably well satisfied 
even for rather thin walls considered in this paper. Therefore we expect that the sources that 
are not captured by the gradient expansion Eq. (Tl9l) — an important example of which is 
quantum mechanical CP-violating reflection [2| — are subdominant. This then suggests that 
the source (Tl9l) provides quite generically the main contribution to the chargino-mediated 
BAU in the nMSSM. 

In the MSSM, where q(z) = and d z \i = 0, the evaluation of the second-order source (fl9"j) 
leads to 

sSU-M ! .^-HJ: (20) 

using the conventions of Ref. [9j , this can be written as 

Tr sSsm ~ 9 2 (0;> 2 + Ml) d k J e ^\ (21) 

where A denotes the difference of the eigenvalues of the matrix rrvm. 

In the nMSSM, there are various additional contributions from the derivatives acting on 
H in the source ( Tl9|) . especially a novel diagonal term of the following form: 

S$ssm~< l „ ° ■ ( 22 ) 

This contribution dominates if fi <C M2, which is usually the case in the nMSSM since 
fi is related to the singlet vev, while the Wino mass parameter M 2 is not related to the 
parameters of the singlet sector and expected to be of the SUSY scale. 
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Hence, we will consider two scenarios. First, we consider the case of large M 2 ~ 1 TeV. 
In this region the contribution (122]) will almost coincide with the full expression (1191) . Sec- 
ond, we will choose a rather small Wino-mass parameter M2 ~ 200 GeV. In this case the 
additional contributions in Eq. ( fl9i) can lead to an enhancement or a cancellation in the 
BAU, and one should keep in mind that the neglected mixing effects could contribute as 
well. In both cases, using the full second-order source (fT9|) . the baryon to entropy ratio r\ is 
determined as was done in Ref. [§] for the MSSM. 

We use a system of diffusion equations that was first derived in Ref. [34| and later adapted 



in Refs. 



6 



7 



11 



35]. This system describes how the CP violation is communicated 



from the chargino sector to the left-handed quarks and finally biases the sphaleron processes. 
These diffusion equations rely on certain assumptions, e.g. that the supergauge interactions 
are in equilibrium, and hence lead to sizable uncertainties. Furthermore, we do not take 
into account recent developments in the determination of interaction rates presented in 



Refs. 



36 



371 ] . but employ the parameters of Ref. p]. Nevertheless, the accuracy of the 



determined BAU should be sufficient for the analysis in this work. 

In the following we will briefly focus on the term ( l22i) . which is prominent in the limit 
of large M 2 . Beside the critical temperature T c , the generated BAU only depends on the 
profile of the Higgsino mass parameter fi(z) during the phase transition, namely the change 
of the phase Aq s , the wall thickness l w and the profile of |/x(z)|. 

With good accuracy, the baryon to entropy ratio, rj = n^/s, scales as 

Ms 



7] OC 



/ T 



(23) 



The dependence on the profile of \\i{z)\ is shown in Fig. [H In this example, the profile is 
parametrized by 



\fi(z)\ = Ho - A/j, ( ~ + - tanh (z/k 



(24) 
(25) 



and the values /j,q = A//, l w = 10/T C , Aq s = n/10 and several values of T c have been chosen. 



To give some feeling about the produced BAU, we note that a good estimate of the 
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50 100 150 200 250 300 350 400 450 500 
Au /GeV 



FIG. 1: The dependence of r]i = 10 10 rj on the chargino mass parameter Afi = fj, . The parameters 
used are l w = 10/T C and Aq s = ir/W. 

predicted r] 10 is in this case given by the formula 

3 

Vio - c(T c ) ^-J- (J±) ' ^ exp(- / , /rT c ), (26) 

7T l w l c \T1 C J Tl c 

with c(T c ) pa 1.6T c /GeV and r « 0.78. This formula characterizes the BAU in the case of 
large M 2 « 1 TeV and /i ~ A/i. 



B. EDM constraints in the nMSSM 



The most severe experimental constraints on CP violation come from measurements of 
the EDM of the electron, gL < 1.6 x 10~ 27 ecm 



381 ] . and neutron, d n < 3.0 x 10 26 ecm 



Already at the one- loop level, contributions of the superpartners give sizable effects in the 
case of CP-violating phases of 0(1). These one-loop diagrams are the same in the nMSSM 
and MSSM, setting \i = —\(S). Minimizing the Higgs potential, we can compute the phase 
of the effective \x parameter. Of course, this phase could be neutralized by introducing a 
compensating phase in the parameter M 2 . Such a tuning would allow us to eliminate the 
one-loop EDMs completely, without much affecting the generated baryon asymmetry, since 
the dominating source is proportional to the change in the phase Ag s and not sensitive to the 
value of q s = arg(/x) in the broken phase. This is not possible in the MSSM, because there 
the produced BAU is, like the electron EDM contribution from the charginos, proportional 
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to the combination Im( / uM 2 ). 

However, here we take a different approach, using only the phase q t in the Higgs potential 
as sole source of CP violation. The one-loop EDMs then induce mass bounds for the first 
and second generation squarks and sleptons; depending on the model parameters these are 
in the range from a few TeV up to 50 TeV [40]. Since the constraints from the neutron EDM 
are usually less stringent than the ones coming from the electron EDMs, we will focus in the 
following on the latter. In single cases we also calculated the Barr-Zee type contributions 
to the neutron EDM 41] . but they barely reach the most recent experimental bounds 39]. 



Besides heavy sfermion masses, there is another possibility to suppress the one-loop elec- 
tron EDM in our model. Notice that the absolute value of the CP-violating phase (q s + q) 
can be smaller in the broken phase than in the symmetric phase. This phase (q s + q) is 
the only CP-odd combination that enters the electron EDMs on the one-loop level in our 
simplified nMSSM model, with CP violation only in the t s parameter. Hence, there is the 
possibility that today's observed sin(g s even though the phase q s greatly changed 

during the phase transition thus producing the BAU. We will analyze this possibility in 
detail in Sec. IVI1 The explicit form of the one-loop electron EDM contributions is given in 
Appendix [B] This possibility entails a certain amount of tuning. 

Additional EDMs can be generated from two-loop chargino or Higgs graphs (see, for 



instance, 



40] and references therein). Notice that the MSSM two-loop chargino contribution 



to the electron EDM 



19 



20] is proportional to tan/? and hence subleading in our model 



that usually predicts tan/? ~ 0(1). Potentially harmful diagrams, including the additional 
CP-odd scalar in the nMSSM, are small as well because of the modest tan/3 and because 
only the Higgs component of the CP-odd scalars couples to the charginos, while the singlet 
component delivers the additional CP violation (for a calculation of these contributions see 



V. ELECTROWEAK PHASE TRANSITION 



One of the parameters entering our baryogenesis analysis is the thickness of the Higgs 
wall profile during the electroweak phase transition Since our CP- violating source is a 
second-order effect in gradients, the integrated BAU scales as 77 ~ l/l w , as already mentioned 
in Eq. ([23]). 
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To determine the wall thickness l w one has to examine the dynamics of the phase tran- 



sition 



This has been done for the MSSM in Ref. [12| and for the NMSSM in Ref. [26|. 
Typical values for the MSSM seem to be close to l w = 10/T C . In the nMSSM, we expect 
rather thin wall profiles, since the linear singlet term and the trilinear singlet Higgs term in 
the effective potential will make the phase transition much stronger than the loop-suppressed 
stop corrections that are responsible for the first-order phase transition in the MSSM. In 
the light of Eq. fl23|) this will further enhance the produced asymmetry with respect to the 
MSSM case. 

To determine the dynamical parameters of the wall, we solved for the classical bounce 
solution of the Higgs and singlet fields ((pi,<p2,<Ps,Q,qs) at the critical temperature (where 
the two minima of the potential are degenerate). Our numerical approach is based on the 
variation of the classical action and is discussed in detail in Ref. j^J]. 

Another parameter relevant to the dominating source in Eq. (122]) is the profile of the CP- 
violating phase of the singlet field. A typical solution is displayed in Fig. [2] corresponding 
to the parameters Aq s = 0.119 and l w = 4.81 T~ x . 

-0.76 | , , , , r- 

-0.78 
-0.8 
-0.82 

cr 

-0.84 
-0.86 
-0.88 

-0.9 ' 1 1 1 1 L 



ArgS 



10 20 30 40 50 60 70 80 90 100 



z/T c " 1 



FIG. 2: A typical wall profile for the parameter q s , corresponding to the parameters Aq s = 0.119 
and l w = 4.81 T" 1 . 



To illuminate a little bit the nature of the phase transition, we will recall some approxi- 
mate analytical criterion for the occurrence of a first-order phase transition, first given in a 
slightly different way in Ref. 18) . Consider the tree-level potential in Eq. ([6]) without CP 
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violation, q t = 0. Additionally assume that the temperature effects give rise to the effective 
potential 



AV T = a^ 2 T 2 , 
where a is some unspecified positive constant. In Ref. 



(27) 



181 ] it was shown that a necessary 



condition for a first-order phase transition is approximately given by 



m. < — 
A 



m R 



m s a 



This can be seen in the following way. Given a certain value for 
evaluated to be 



(28) 

can easily be 



t s + a 



m? s + A 2 2 

Using this in our potential, and expanding around 



t 



V + AV T = -^+(M 2 + aT 2 )0 2 + A 2 4 - 



0, we obtain 

(t a + ~a<p 2 ) 2 



C + Ci(j) 2 + C 2 4 + C 3 6 



m 2 + \ 2 (f) 2 



with the coefficients 



Co 



Cl 



c 2 



c 3 



tz 



mi 



M + aT 
1 



2at« 



mi 



+ 



2+2 

4 ; 



A 2 t 



A - 

A^ 

,4 



ml 



A 2 t. 



a — 



mi 



\ 2 U 



(29) 



(30) 



(31) 



mi \ m% 

If the symmetric minimum is absent at zero temperature, Ci(T = 0) < 0, a temperature T 2 
can be found, such that ci(T 2 ) = 0. For this temperature there exists a lower-lying potential 
minimum in the case c 2 < 0, which is equivalent to the condition in Eq. ([25]) . Since for 
temperatures T > T 2 a potential well develops between the symmetric and the lower broken 
vacuum, a first-order phase transition is possible, given the vacuum decay rate is large 
enough such that the transition occurs before the temperature T 2 is reached. Hence, it is 
possible in the nMSSM to obtain a first-order phase transition due to tree-level dynamics, 



16 



in contrast to the MSSM. Analyzing the numerical results, we will see that the constraint 
in Eq. (1251) is usually fulfilled in viable models, even if the one-loop contributions to the 
potential and the CP phase are included and hence that the phase transition is dominated 
by the tree-level dynamics. 

The argument just presented is a concrete realization of the general effective field theory 



approach recently discussed in Refs. 45j, |46j. There it was shown that a strong first-order 
phase transition can be induced at tree-level by the interplay of a negative 4 term and a 
positive 6 term which stabilizes the potential. The suppression scale of the 6 term should 
be somewhat below a TeV for the mechanism to work. Here the relevant 4> 4 and 6 operators 
are generated by integrating out the singlet field. This generalizes the usual situation, where 
the phase transition is induced by a negative 3 term and a positive 4 term. 

VI. NUMERICAL ANALYSIS 

To inspect the parameter space of the nMSSM we proceed as follows. First, we choose 
random parameters for the Higgs potential in the ranges displayed in Tab. [TTJ To ensure 
maximal numerical stability, all chosen parameters are of 0(1), and can be thought of as 
dimensionless parameters. Those parameters then lead not to the physical Higgs vev, but 
to some dimensionless Higgs vev 0o- Finally, all dimensionful quantities, such as the critical 
temperature or the mass spectrum have to be scaled with (173.458 GeV)/0o to yield the 
physical values. During the minimization of the potential, stable and metastable broken 
minima were analyzed. Depending on the parameters, metastable minima occur but, in 
the CP-conserving case q f = 0, no transitional (spontaneous) CP violation was observed in 



contrast to the NMSSM |26|] that contains an additional cubic singlet term. 

Next, we correct the bare parameters with the one-loop contributions of Eq. (jHJ and 
confront them with the constraints on the mass spectrum from Tab. [T] and on the Z-width 
from Eq. (1T61) . 

If the parameter set passes these constraints, we add the temperature dependent contri- 
butions to the effective potential as explained in Sec. Ill Cl and examine the phase transition. 
We require that the models have a first-order phase transition of sufficient strength H] , 
<Po/T c > 0.9. 

Before we discuss baryogenesis in our model, we would like to examine restrictions on the 
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lower bound 




parameter 




upper bound 


0.1 


< 


m\ 


< 


1 


0.1 


< 


m 2 


< 


1 


0.1 


< 


mi2 


< 


2 





< 


m s 


< 


2 


-2 


< 


X 


< 


2 





< 


OA 


< 


2 





< 


1/3 


< 


2 





< 


qt 


< 


2vr 



TABLE II: Dimensionless parameter ranges used for the numerical analysis. 



parameters imposed by the constraints on the mass spectrum and comment on the criterion 
for a first-order phase transition given in the last section in Eq. (I28I) . 

First, the eight parameters given in Tab. [II] have to lead to the correct Higgs vev, which 
is achieved by a rescaling of the dimensionful parameters. Hence our parameter space is 
effectively only seven dimensional. One restriction on the parameters is that a\ cannot be 
chosen arbitrarily large, since this destabilizes the potential in the negative S direction. 
Analyzing the parameter sets that fulfill the mass constraints, one observes that the pa- 
rameters A and q t are not distributed homogeneously. Small values of A make it seemingly 
difficult to fulfill the mass bound of the chargino, since one of the diagonal entries of the 
chargino mass matrix is — A0 S . To have a potential with extremely large vev <p s requires at 
least some fine-tuning, since the one-loop contribution tends to yield an effective potential 
that is unbounded from below, if the tree-level parameters are chosen to provide a large vev 



4> s - Large values o: 
to a Landau pole 



A hence seem to be the more natural choice, even though they can lead 



isj]. Usually, the mass constraints on the neutralinos are automatically 
fulfilled, if the charginos surpass their more restrictive bounds, but additional constraints on 
the parameters enter through the spectrum of the Higgs particles. In many cases, a range 
of values for the parameter a\ can be found, where off-diagonal elements in the Higgs mass 
matrix cancel, which tends to enlarge the lightest Higgs mass. In addition the parameter a\ 
has a strong influence on the phase transition according to Eq. (128]) . 

This situation is demonstrated in Fig. [3]for a parameter set with a rather small parameter 



18 



160 I 1 1 1 1 1 1 1 r- 

150 [ .. T c 

brok 



Vi 



140 
130 
| 120 
110 
100 
90 
80 



172 173 174 175 176 177 178 179 180 181 182 

/ GeV 

FIG. 3: The critical temperature T c , the Higgs vev <f> in the broken phase at T c and one Higgs mass 
as functions of a\. 

A. The parameters mi, m 2 , m 12 , m s are chosen such that tan(/3) = 2.0, S = —250 GeV, 
= 173 GeV and M a = 500 GeV at the tree- level, where the CP-odd Higgs mass parameter 
is defined by 

M 2 a = . 1 {m\ 2 - a A s ) . (32) 
sin p cos p 

The remaining parameters are A = 0.55, t] /3 = 70 GeV, q s = 0.3, while a\ is varied. For 172 
GeV < a\ < 178 GeV, this model develops a strong first-order phase transition and generates 
more than the observed BAU. For lower values of a\, the model has no stable broken phase, 
while for larger values of a\, the phase transition is too weak. The plotted Higgs mass is 
that of the third lightest Higgs, but the two lighter states would have escaped detection at 
LEP because of the suppressed coupling to the Z-boson. This example demonstrates that 
even though there exist viable models without Landau pole in A, this possibility entails a 
certain amount of tuning. 

The second parameter that is restricted by the mass constraints is the CP- violating phase 
q t . The reason for this effect is that values with cos(q t ) ~ —1 lead to smaller Higgs masses. 
Fig. H] displays the parameters A and qt for a set of random models that fulfill the mass 
constraints. 

Demanding a strong first-order phase transition further restricts the parameter space. In 
Fig. [5] we plot the left-hand side versus the right-hand side of the criterion in Eq. ([28]) (both 
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FIG. 4: The parameters A and q t for a set of random models that fulfill the mass constraints. 

sides are scaled by l/(m 1 m 2 ) to make them dimensionless) . In the left plot we use random 
models that fulfill the mass constraints on the spectrum, but are unconstrained otherwise; 
in the right plot we impose the mass constraints on the model and require a strong first- 
order phase transition. In the latter case, most of the parameter sets are in accordance 
with Eq. ( l28l) . while the parameter sets in the former case are evenly distributed. Hence, 
the tree-level criterion for the phase transition seems to be applicable even if the one-loop 
contributions to the effective potential and the CP phase are taken into account. 

In the nMSSM, for several reasons, we expect a much larger BAU than in the MSSM. 
First, the parameter tan(/3), which needs to be large in the MSSM, is naturally 0(1) in the 
nMSSM as depicted in Fig. El This does not only help to suppress the two-loop contributions 
to the EDMs, as discussed in the previous section, but also enhances the contributions from 
the source in Eq. (|21~1) . Secondly, the wall thickness, which in the MSSM usually is of order 
20/T c -30/T c |l2j, can be much smaller, since the phase transition is strengthened by the 
linear and trilinear terms in the effective potential. This claim is supported by Fig. [7J The 
third reason for the enhancement of the BAU in the nMSSM with respect to the MSSM is 
the additional source in Eq. ( 1221) . which is in many cases dominating the generation of the 
BAU. 

Finally, we calculate the generated baryon asymmetry and compare the result with the 
experimental observation, rj = (0.87 ± 0.03) x 10 -10 48j. The result is shown in Figs. [8] 
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and [9j Approximately 50% of the parameter sets predict a higher value than the observed 
baryon asymmetry in the model with large M 2 = 1 TeV, while this number increases to 63% 
for small M 2 = 200 GeV. 

In addition, we plotted the BAU generated by parameter sets that fulfill the experimental 
bounds on the electron EDM with sfermion masses of 1 TeV in the first and second gen- 
eration. Some of them predict a BAU in accordance with observation, and hence give the 
possibility to construct nMSSM models that contain less constrained sfermions (lighter than 
1 TeV) being at the same time consistent with EDM constraints and baryogenesis. In some 
cases the electron EDM is small because of a random cancellation between the neutralino 
and chargino contributions, but occasionally the suppression of the electron EDM is due to 
the fact that the combination sin(g s + q) is relatively small in the broken phase. 



VII. CONCLUSION 

We have analyzed the phase transition and baryogenesis in the nMSSM (TjQ) - (j4j) with CP 
violation in the singlet sector. We have shown that the singlet field enhances the strength 
of the phase transition in such a way that one typically obtains a strong phase transition, 
as required for successful baryogenesis. This is to be contrasted with the MSSM, in which 
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FIG. 5: The plots show the combinations of the parameters that enter the tree-level condition 
for a first-order phase transition, Eq. (|28p . The left plot contains parameter sets that fulfill only 
the mass constraints, while the right plot contains parameter sets that have in addition a strong 
first-order phase transition. 
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FIG. 6: The binned result of the parameter tan(/3) for the parameter sets with a strong first-order 
phase transition. 



the mass of the lightest Higgs field must not be greater than about 120 GeV, and the 
right-handed stop must be light, mf G (120, 160) GeV. 

Next we performed the calculation of baryogenesis mediated by charginos in the nMSSM. 
After calculating the CP- violating sources in the gradient expansion, we argued that in most 
of the parameter space the dominant source comes from the second-order semiclassical force 
in the Boltzmann transport equation for charginos. The source related to flavor mixing, of 
first order in the gradient expansion, tends to be smaller, because the bubble wall is rather 
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FIG. 7: The binned result of the parameter l w for the parameter sets with a strong first-order 
phase transition. 
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FIG. 8: The binned result for the BAU analysis for large Mi = 1 TeV. Approximately 50% of 
the parameter sets predict a value of the baryon asymmetry higher than the observed one. The 
bottom line corresponds to parameter sets that fulfill current bounds at the electron EDM with 1 
TeV sfermions (4.8%). 

thin. For a generic choice of parameters, one is far from the chargino mass degeneracy, where 
the first-order source may be important. To come to this conclusion, we used an approach 
to the calculation of the first-order sources 

aa 

that differs from earlier work 
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FIG. 9: The binned result for the BAU analysis for small M2 = 200 GeV. Approximately 63% of 
the parameter sets predict a value of the baryon asymmetry higher than the observed one. The 
bottom line corresponds to models that fulfill current bounds at the electron EDM with 1 TeV 
sfermions (6.2%). 
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in the sense that our treatment of sources is basis independent, and the magnitude of the 
source in the transport equation is unambiguous. Using this advanced transport theory, 
the first-order sources are of a somewhat lower amplitude and exhibit a much narrower 
resonance near the chargino mass degeneracy, with the effect that in most of the parameter 
space the second-order source dominates. In the MSSM, this is not the case for the chargino- 
mediated baryogenesis because the wall tends to be thicker, thus weakening the second-order 
(semiclassical force) source, while leaving the first-order source more or less unchanged. 
Furthermore, the dominant second-order source of baryogenesis in the nMSSM, Eq. (1221) . is 
not present in the MSSM. Owing to these differences, successful baryogenesis in the MSSM is 
only possible near the resonance (chargino mass degeneracy), and with nearly maximum CP 
violation, which is in conflict with the current EDM bounds, unless a tuning of parameters 
is invoked jc|. 

On the other hand our analysis of the baryon production in the nMSSM looks promising. 
When we restrict the CP-violating phase in the singlet sector to be about qt ~ 0.3, and take 
tan(/3) ~ 1, we still get baryon production consistent with the observed value. This choice 
certainly does not violate any of the current EDM bounds. When q t is chosen randomly, 
approximately 50% (63%) models predict more than the observed BAU when M 2 = 1 TeV 
(M 2 = 200 GeV), which indicates that baryogenesis in the nMSSM is generic. 

It is finally interesting to compare baryogenesis in the nMSSM and in the general NMSSM 



formerly analyzed in Ref. 26]. Because of the presence of a singlet self coupling and an 
explicit \x term, the NMSSM allows for a much richer Higgs phenomenology. (However, 
some additional assumptions about the structure of the higher dimensional operators have 
to be made, in order to prevent the destabilization of the electroweak scale by corrections 
to the singlet tadpole.) There is no danger of a light singlino state, so that the model can 
account for the observed baryon asymmetry also for values of A < 0.1 and large tan/3. It 
also allows for transitional CP violation, which means an electroweak phase transition that 
connects a high-temperature CP-broken phase with a low-temperature CP-symmetric phase. 
This way there are varying complex phases in the bubble wall, without leaving a trace at low 
temperatures. In the present model this possibility is prevented by the strongly constrained 
Higgs potential. 



In Ref. 



26j the source terms were computed by solving the Dirac equation for the 



charginos in the WKB approximation. As has been recently shown in Ref. H[ , this formal- 
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■ B, 



ism can reproduce the second-order source of Ref. [5|, used in the present work, if the Lorentz 
transformation to a general Lorentz frame is done carefully (see also Ref. 49J). These results 
suggest that in Ref. [261] the baryon asymmetry was underestimated by a factor of 2 to 5. It 
seems to be interesting to update this analysis to meet the present experimental constraints. 
It also seems promising to apply the presented techniques to more general supersymmetric 
models, such as models with extra U(l) symmetries 501 ]. 
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APPENDIX A: MASS SPECTRUM OF THE NMSSM 

In the following we collect the mass matrices that have been used in the one-loop poten- 
tials in Sec. flU 

We used the physical constants 

sia(9 w ) = 0.2312, a EW = 1/127.907, (Al) 
m z = 91.1876 GeV, m t = 165.0 GeV, (A2) 

which give rise to the values 

g = 0.357, g' = 0.652, O = 173.458 GeV, (A3) 

where O denotes the T = value of the vev 0. 

If not stated differently, we have used the SUSY-breaking parameters 

mu = mQ = 500 GeV, a t = 100 GeV, 

M 2 = 2M X = 200 GeV, m E = 500 GeV, (A4) 



which enter in the mass matrices that we will define subsequently. 



1. Higgs bosons 



For the neutral Higgs bosons we use the notation 

H? = e^ifa + S t + iPi), S = e^(0 s + S s + iP s ), 
and q\ = q 2 = q/2. The corresponding mass matrix has the following form 



2 | Mis M 2 SP 



M PS m pp 



with the matrices M S s, M S p = M PS , and M PP given by the following entries. The CP- 
entries at the tree-level read 

Mlsn = ro? + A 2 (0 2 + + ^(30 2 - <f>l), 

M lsi2 = ~ m i2 cos (<?) + °A0 S cos(g + q s ) + ^0i0 2 (4A 2 - g 2 ), 
M ssi3 = °A02 cos(g + q s ) + 2A 2 0i0s, 

M 2 SS22 = ml + \\4>\ + + 9 ^4>l ~ 4>\), 
M SS23 = a A 0!Cos(g + g s ) +2A 2 2 S , 

M SS33 = m l + A 2 2 - 

The CP-odd entries are 



7,2 



M 2 PPU = m\ + A 2 (0 2 + 2 ) + ^(0 2 - 2 ), 

Mf, P12 = m\ 2 cos(g) - a A s cos(g + q a ), 

m pp\3 = -a A 02COs(g + q a ), 

q2 

M 2 PP22 = m\ + A 2 (0 2 + 2 ) + ^(0 2 - 2 ), 

m pp23 = -a A 0icos(g + g s ), 

2 _2 , \2/2 



m PP33 = mi + y<p z . 

Finally the CP-mixed entries yield 



M SP12 = m 2 u sm(q) - a X (p s sm(q + q s ), 
M SP13 = -a X (f) 2 sm(q + q s ), 
M SP23 = -OA0isin(g + q s ). 
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In the CP-conserving case the submatrix (1A9I) vanishes so that CP-even and CP-odd states 
do not mix. 

If the one-loop effective potential is included, we determine the masses of the neutral 
Higgses by the second derivatives of the effective potential. 

The mass matrix of the charged Higgs bosons in the basis (H^ , H£) contains the complex 
entries 

^2 „2 



n \\ 4 A 



M h = -2^2(2A 2 -/) + m^-a A s e 

^2 



i(q+q a ) 



M 



ml + X^- 9 -( ( pl-4>l) + 9 - ( j ) l 



(A10) 



2. Charginos and Neutralinos 



The chargino mass matrix reads 



yL 







Mo 



!> 2 e 1 2 



gcfiie *2 — \(j) s e iqa 



and is diagonalized by the biunitary transformation Y. 



diag 



U*Y*±V*. 



(All) 



The symmetric mass matrix of the neutralinos in the basis x° = (iB,iW 3 , H®, H®, S) 
yields 

. A 



M x° 



Mi 









M 2 















\<p s e iq ° 








A0 2 e*2 



and can be diagonalized using a unitary matrix, M~^ s = X T M^oX. 



(A12) 



3. Gauge bosons 

The mass of the PF-boson is given by 



m 



^ 2/2 



(A13) 



27 



while the photon and the Z-boson share the following hermitian mass matrix 

Ml = \ 2y v 2 yy ^ A14) 

that leads to the Z-boson mass mz = (g/y/2)(j). 

4. Tops and stops 

The top quark has the mass 

m 2 t = yUl (A15) 
and the masses of the stops are given by the following hermitian matrix 

f m 2 Q + m 2 + \{g 2 - \g' 2 ){<P\ - <j%) a t <f) 2 e^ 2 + y^fae^^ 2 ' 



Mf=\ v 1 4W 6 ° /vrx rz/ — (A16) 

1 atfce-*' 2 + y t X(j)sM i(qs+q/2) m 2 + m 2 + \g' 2 {^\ - <f> 2 2 ) 



5. Sneutrinos and selectrons 

For the selectrons we have 

/ rriE —Xd> s e iqB m e tan(B) \ 

Mj=\ . £ e W , (A17) 

\ — A0 s e i9s m e tan(/?) m E J 

that is diagonalized via the transformation M| diag = M 2 D. 

For the electron EDM contributions from the charginos we use a sneutrino of mass = 
1 TeV. 



APPENDIX B: ONE-LOOP CONTRIBUTIONS TO THE ELECTRON EDM 

In this Appendix we briefly discuss the one-loop contributions to the electron EDM 
coming from chargino and neutralino exchange. For the sfermions of the first twogenerations 



we assume masses of 1 TeV. The contribution from the charginos is given by [511 ]: 

4-cha rK ino a EM yl^j^.) A r^Y (B1) 



i=l 

where 



4irsm 2 (6 w ) _ 

iq/2 



T ei = K e U* 2 V a , n e = me€ — . (B2) 

\/2mw cos(p) 
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Here, the matrices U and V diagonalize the chargino mass matrix as defined in the previous 
section and 



A(r) 



1 



3 — r 



2(1 -r) 2 V" ' 1 - r 

denotes the loop function. 

Analogously, the contribution from the neutralinos is 

2 4 



ln(r* 



(B3) 



51 



52] 



d, 



e— neutralino 



/e 



Air sin 



olem 

27" 



7 w 



k=l i=l 



EE^M^) B(m|o/ 



(B4) 



with 



(tan{d w )X li + X li )D* k + V2K e X 3i D* k 



tan{6 w )X u D* 2k + — X u D\ k 



(B5) 



In this equation, X and D diagonalize the neutralinos and selectrons, respectively. The loop 
function B is defined by 

1 / 2r 



Br 



Following the discussion in Ref. 



2(1 -r) 



1 +r + 



1 — r 



ln(r' 



(B6) 



51| , it can be shown that if q t is the only CP- violating phase 



in the Higgs sector, both contributions depend only on the CP-odd combination 
The current experimental bound on the electron EDM is d e < 1.6 x 10 _27 ecm 



APPENDIX C: EXAMPLE SETS 



In this appendix we give some examples of parameters that develop a strong first-order 
phase transition. The examples are not chosen arbitrarily, but they represent specific cases. 
The first set, shown in Tables [IIJHXl accomplishes the generation of a large BAU thanks to a 
relatively thin wall. On the other hand the electron EDM is extremely small, partly because 
the combination q + q s is rather small in the broken phase and partly due to a coincidental 
cancellation between the neutralino and the chargino contributions to the electron EDM. 
Set number 6 describes a model with small A, taken from Fig. [3] with a\ = 177.3 GeV. This 
model generates more than the observed BAU, but the calculated electron EDM is slightly 
too big, so that the sfermions of the first two generations have to be heavier than 1 TeV. 
Starting from q t ~ 0.1 instead of 0.3 should yield r/io ~ 1 and an electron EDM within the 
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experimental bound. The remaining sets are randomly chosen. Set number 6 demonstrates 
that even a sizable value of q t does not lead to a large baryon asymmetry if Aq s is small 
along the bubble wall. The rather thick wall induces an additional suppression of 7710 in this 
case. 
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set 


mi in GeV 


m 2 in GeV 


m\2 in GeV 


m s in GeV 


A 


a\ in GeV 


t 1 / 3 in GeV 


qt 


1 


157.2 


93.0 


170.6 


55.8 


1.2642 


268.2 


98.4 


2.113 


2 


124.2 


149.2 


215.0 


127.2 


1.4320 


254.1 


152.2 


5.050 


3 


248.9 


230.8 


243.3 


160.5 


-0.8937 


214.2 


190.5 


0.046 


4 


397.1 


251.9 


375.5 


342.1 


-1.2991 


292.0 


310.5 


0.282 


5 


66.4 


98.5 


111.8 


65.5 


0.7656 


95.3 


102.0 


6.193 


6 


425.1 


165.7 


240.1 


26.8 


0.5500 


177.3 


70.9 


0.300 



TABLE III: Parameter examples used for the numerical analysis. 



set 


X ± 


H± 


h/2 


1 


221.77 


107.40 


219.56 


521.43 


529.30 


2 


221.64 


131.22 


236.41 


537.16 


514.91 


3 


270.45 


109.31 


410.84 


485.74 


563.03 


4 


146.21 


315.37 


588.45 


479.29 


567.75 


5 


221.06 


105.94 


148.97 


530.64 


521.67 


6 


225.24 


144.62 


494.87 


520.50 


528.89 



TABLE IV: The mass spectrum of the charginos, charged Higgses and stops (in GeV). 



set 


S,P 


1 


142.72 


210.53 


217.41 


273.78 


357.83 


2 


177.85 


260.59 


292.50 


296.78 


396.55 


3 


119.12 


202.13 


251.16 


432.44 


455.77 


4 


121.65 


395.09 


448.97 


602.98 


641.00 


5 


115.07 


118.67 


169.87 


211.08 


237.02 


6 


76.74 


89.86 


114.76 


504.68 


506.91 



TABLE 



V: 



The mass spectrum of the neutral Higgses (in GeV). 



set 


X° 


1 


267.76 


105.80 


113.09 


221.89 


184.03 


2 


313.48 


105.14 


138.79 


222.79 


198.12 


3 


135.40 


68.73 


89.48 


275.43 


269.04 


4 


389.60 


323.06 


81.50 


161.87 


123.56 


5 


222.16 


181.83 


116.61 


107.33 


94.72 


6 


227.33 


42.25 


105.37 


165.39 


175.04 



TABLE VI: The mass spectrum of the neutralinos (in GeV). 



set 


4> in GeV 


P 


Q 


<p s in GeV 


Qs 


1 


173.46 


0.926 


0.157 


-71.0 


-0.392 


2 


173.46 


0.739 


-0.214 


-81.8 


0.687 


3 


173.46 


0.808 


0.015 


-202.9 


-0.036 


4 


173.46 


0.915 


0.074 


-202.6 


-0.256 


5 


173.46 


0.714 


-0.024 


-112.6 


0.051 


6 


173.46 


1.108 


0.029 


-251.9 


-0.067 



TABLE VII: The vevs in the broken phase at temperature T = 0. 



set 


(f> in GeV 


(3 


<? 


4> s in GeV 


<ls 


1 


165.4 


0.900 


0.175 


-71.4 


-0.437 


2 


170.0 


0.726 


-0.223 


-83.2 


0.708 


3 


150.3 


0.803 


0.016 


-214.9 


-0.038 


4 


164.5 


0.911 


0.076 


-207.4 


-0.259 


5 


151.6 


0.687 


-0.029 


-120.2 


0.058 


6 


141.9 


1.102 


0.041 


-261.5 


-0.092 



TABLE VIII: The vevs in the broken phase at temperature T = T C . 



set 


(j) in GeV 


& 


q 


4> s in GeV 




1 


0.0 






-281.4 


-2.113 


2 


0.0 






-241.6 


1.233 


3 


0.0 






-283.3 


-0.046 


4 


0.0 






-259.9 


-0.282 


5 


0.0 






-266.0 


0.090 


6 


0.0 






-471.5 


-0.300 



TABLE IX: The vevs in the symmetric phase at temperature T = T c . 



set 


T c in GeV 


d e in 10 27 e cm 


l w in T~ l 


Vw 


1 


113.5 


0.002 


2.43 


29.443 


2 


99.1 


0.796 


6.38 


-3.214 


3 


109.1 


0.499 


7.82 


0.014 


4 


78.6 


5.894 


34.01 


0.005 


5 


115.8 


-0.893 


3.05 


-0.398 


6 


105.6 


-2.054 


2.39 


2.717 



TABLE X: The parameters of the phase transition and the generated BAU. 
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